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Abstract. We consider a two-dimensional bi-layered loop model with a certain 
interlayer coupling and study its spectrum on a torus. Each layer consists of an 0(n) 
model on a honeycomb lattice with periodic boundary conditions; these layers are 
stacked such that the links of the lattice intersect each other. A complex Boltzmann 
weight A with unit modulus is assigned to each intersection of two loops each from 
■ each layer. The model is reduced to an inhomogeneous vertex model at a special 

point of parameters. The continuum partition function is represented, based on the 
idea of the Coulomb gas, by a path integral over two compact bosonic fields. The 
modular invariance of the partition function follows naturally. Further, because of 
O | the topological nature of the interlayer coupling, the fluctuation of loops decomposes 

into a local and a global part. The existence of the latter leads to a sum over all 
the pairs of torus knots, which can be Poisson ressummed by the Mobius inversion 
formula. This reveals the operator content of the theory. The multiplicity of each 
operator is explicitly given by a combination of two Ramanujan sums. We calculate 
each scaling dimension as a function of A. We present the flow of dimensions which 
' connects the decoupled-O(l) models at A = 1 and the layered-O(l) model with the 

non-trivial coupling A = — 1. The lower spectrum in the latter model is related to that 
of a known coset model. 

O 
O 

> : 

•i-h . 1. Introduction 

X 

H 

Statistical models in two dimensions can often be dealt with by several of the exact 
methods and have been invaluable sources of ideas in quantum field theory and many- 
body physics. However, little is known about two-dimensional models with layered 
structure, which have been considered in various contexts [U [2j [3j HJ E] . In this paper, 
we propose a bi-layered 0(n) loop model with a particular type of the interlayer coupling 
which allows the exact calculation of the partition function by a simple path-integral 
based on the Gaussian action. 

Our main tool is Coulomb gas (CG) method. Let us recall some general features of 
CG. As is well known, fluctuations of a system become important in the lower dimensions 
and, the critical exponents differ much from those of the Gaussian (or mean-field) theory. 
Indeed, in two dimensions, we have a zoo of critical models which are quite different from 
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the Gaussian; for instance, they are the unitary minimal models [TJ, which satisfy highly- 
nontrivial bootstrap conditions under the infinite-dimensional conformal symmetry [6]. 
Interestingly, however, it turned out that introducing intrinsic charges in a Gaussian 
theory can give heuristic, nevertheless exact, descriptions of a certain class of critical 
models pHl CHI 12] • Such descriptions are generically referred as the Coulomb gas and 
enable us to calculate the scaling dimensions [10J, iV-point functions [TTj . and partition 
functions [T2] . 

Let us turn to more specific lattice models whose continuum limits are described 
by the CG. The simplest is the 0(n) loop model (see [TH [15] for recent reviews). This 
model actually have considerable importance for general models with the central charge 
c ^ 1. For instance, each member of the unitary minimal models can be represented as 
a superposition of several different loop models 0- 

In a standard, algebraic view, partition functions of conformal field theories in 
general are assumed to have the sesquilinear form 



where h is the allowed highest weight and Xh is the character of the representation in 
some algebra [jj; an integer N h ^ is the multiplicity of the corresponding representation. 
This expression ([1]) reveals the operator content, from which one can extract various 
physics in the continuum limit. 

The CG method allows us to evaluate the partition function of the 0(n) loop model 
on a torus in a remarkably direct way (without algebraic inputs). It is essentially given 
by a path-integral over a compact bosonic field governed by the Gaussian action 



where the bosonic field (p is a continuum version of the height function on the dual 
lattice which is uniquely determined from the non-intersecting loop configurations on a 
honeycomb lattice [12J . The advantage of the CG methods is this kind of flexibility in 
scales compared to the other exact methods; from microscopic (lattice) formulation of 
a model, we can derive the corresponding macroscopic (continuum) behavior, which is 
yet straightforward to handle analytically. 

In this paper, we extend the 0(n) loop model to a bi-layer model through the lattice 
formulation and calculate its continuum partition function on a torus in the form by 
which one can, in general, read off the operator content. The model has an interaction 

X More concretely, along the lines of [I7J US US] , one can see the partition functions of the 0{n) model 
which is generalized to have another parameter h (a weight for non-contractible loops; see Section [3. ip 
form a complete basis for those of the A-D-E lattice models (| ^ c ^ 1) 17 and therefore for the 
minimal conformal field theories [71 1201 [51]. Further, the CG analysis of the loop models also provide us 
with the results on the non- unitary models such as the percolation, self-avoiding walk (c = 0) [T2l I22j . 
and dense polymer (c = —2) [T2l |2"5] . 

§ The character formulas, for example, relevant to the unitary minimal models is that of the 
representation of the Virasoro algebra [TS] . 




(1) 
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determined from the number of intersections between two loops each from each layer. 
The decoupling of the local fluctuations from the topological ones, which is essential in 
the single-layer loop model p2], is taken over to this model. This is actually due to the 
topological nature of the inter-layer interaction. 

As we mentioned, multi-layered models appear in many places. In particular, it is 
important for the disordered systems pfl [2], [3]. It is the essence of the replica method 
that taking M — > limit gives the physics of disordered system, where M is the number 
of layers. Although it should be noted that our interlayer coupling might differ from 
that arise from the replica method, we think this as a step towards the understanding 
of the replica limit. 

As the principle of the path-integral dictates, we should sum over all the topological 
sectors. Accordingly, we obtain a sum over all the possible pairs of torus knots; each 
knot is characterized by two coprime integers and corresponds to a homotopical class 
of non-contractible loops on the one of the torus-layers. From this expression, we can 
easily see the modular invariance [8] of the partition function. 

To deal with the sum over the coprime integers, we apply the Mobius inversion 
formula [39]. Our expression for the multiplicity in the bi-layer system involves the 
product of the Ramanujan sums, and looks similar to the one which occurs in the 
additive number theory . It contains highly non-trivial selection rules. To illustrate 
this, we analyze the partition function of the bi-layered 0(1) models. 

The rest of the paper is organized as follows. In Section [27TI the model is defined on 
the lattice with doubly-periodic boundary conditions. At a special point of parameters, 
the model becomes an inhomogeneous vertex model. We derive, in Section 12.21 a 
determinant formula which counts the number of intersections. In Section [3TTT we review 
the ideas used in calculating the partition function of the 0(n) model. In doing so, we fix 
the notations of the parameters used in the continuum formulation. Section [3T21 contains 
our main result. We first show the modular invariance of the partition function. Next, 
we evaluate it by the repeated applications of the Mobius inversion formula. The Poisson 
ressumation can then be performed, resulting in an expression for conformal dimensions 
of various operators, which depends continually on the interaction parameter T. Section 
14. II is devoted to a brief inspection of the examples in the bi-layered 0(1) models. The 
lower spectrum seems to be related to that of a particular coset model. In Section 14.21 
we observe generic properties of the flow of the conformal dimensions when one varies 
T. In Section [5j we conclude with few remarks. In Appendix A , a short derivation of 
the formula concerning the Mobius inversion is given. In Appendix B the critical Ising 
parition function and its squares, which we need for checking one of our examples are 
reviewed from a unified, orbifold point of view. 

|| More specifically, it resembles the "singular series" [JI] (see [351 HE| f° r physical application) which is 
introduced to study both the (asymptotic) Waring's problem and pair-correlation of two prime numbers. 
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2. Lattice partition function of the layered loop model 

2.1. Interlay er coupling depending on the number of intersections 

The 0(n) model is a spin model generalizing the Ising model, whose partition function, 
on a honeycomb lattice, allows the following loop representation: 

Zo(n); lattice = X L n° , (3) 

loops 

where the sum is over non-self-intersecting loop configurations; L and C is the length 
of loops and number of loops, respectively. For \n\ < 2, the model becomes critical at 
some x = x c (n). The universality class of the critical 0(n) model is studied well by 
mapping the model to a six-vertex model on a Kagome lattice with the Coulomb gas 
method [10] or by mapping it to a seven-vertex model on the same honeycomb lattice, 
which is solved by the Bethe ansatz method [37J, EE]- 

A crucial step existing in both methods is to give an orientation to each loop. 
Namely, one decomposes the real weight n in ((3]) into two complex weights: 

n = 2cos7rx = n + + n_, with n± = e ±l7rx , (4) 

and assigns n + (resp. nJ) to each clockwise (resp. anti-clockwise) loop. The partition 
function ((3]) is then written as 

Zo(n); lattice = ^ X L n+ + n^ _ , (5) 

oriented 
loops 

where the sum is over the configurations of oriented loops and C + (resp. C_) is the 
number of clockwise (resp. anti-clockwise) loops. Using the relation |#left turns — 
fright turns | = ±6, which holds for the loops on the honeycomb lattice on a plane, the 
weight n± can be further decomposed and redistributed as: 

e +i7rX//6 (resp. e - " ^ 6 ) on each left (resp. right) turn, (6) 

at a vertex, resulting in the seven- vertex model [371 EE] • 

We define the layered model via the summation ([5]) over the oriented loops. 
Consider two layers A and B each of which consists of the honeycomb lattice with 
doubly-periodic boundary conditions. To obtain an oriented loop configuration, place 
arrows on some of the edges of A and B independently such that at each vertex there 
are either no arrows, or just one arrow in and just one out. Now define a lattice formed 
by stacking A and B such that there are two types of the vertices, V 3 with three edges 
and Vq with six edges; we denote this by Af)B (Fig. [1]). Then two oriented loops I a on A 
and lgon5 can intersect each other on ADB at the vertex of type Vq. As an interlayer 
coupling, we assign a weight A°", (a = ±1) for each intersection of loops depending on 
the orientations. The sign o~p at an intersection P is defined locally and relatively for 
two oriented loops I a and Z#. Namely P is counted as o> = +1 (resp. o P = —1) if l B 
crosses I a at P from right to left (resp. from left to right). 

Let E be the total number of signed intersections: E = Ylii a P l i where Pi is the 
i-th intersection. One easily verifies that E — if there are only loops each of which 
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Figure 1. (a) The stacked lattice AHB formed by the honeycomb lattice of the layer 
A (solid lines) and that of the layer B (dotted lines), (b) A sample loop configuration 
on AnB endowed with the doubly-periodic boundary conditions. This contains one 
non-contractible loop on each layer, and three contractible loops; two on A and one on 
B. The total number of intersections E equals —1. 



is homotopic to a point. Therefore, in order to obtain a coupled model, we have to 
consider the model on a manifold on which some of the loops can be non-contractible 
so that we may have E ^ 0; we here consider the model on a torus. 




Figure 2. Sample of the weights for arrow-configurations at each vertex of the type 
Vq. The solid lines and dotted lines (both with arrows) represent the loop segments 
on the layer A and on the layer B, respectively. These weights are invariant under 
27r/3- rotations of each configuration. Note that the weights at the top and bottom in 
the last column also appear in the set of weights at the vertex of the type V3. 
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(7) 



oriented 
loops 



where L, C + (resp. C_) stand for a total length of loops, a number of contractible 
clockwise (resp. anti-clockwise) loops as before. Since the interlayer coupling is defined 
locally, the model ([7]) can be thought of as a vertex model on the lattice ADB. The 
weight for an oriented-loop configuration on the lattice ADB is given by the product of 
the local vertex weights (Fig. |2]). This vertex model is special in its inhomogeneity: the 
vertex of type V3 allows 7 arrow-configurations and the vertex of type Vq allows 7 2 = 49 
configurations. It would be nice if one could find the exact solution of this and compare 
that with an analysis based on our continuum partition function. 

For an oriented non-contractible loop, we have another relation |#left turns — 
fright turns| = 0, which is distinct from the relation for a contractible loop. Therefore, 
under the vertex rule (EJ), the non-contractible (oriented) loop in (JTj) has the trivial 
weight 1. Instead, we can assign a different weight h = h + + h_ with n + = n*_\ 



where D + and D_ are number of positively and negatively oriented non-contractible 
loops, respectively. For definiteness we use a convention in which we call "positive" for 
a loop of the winding numbers (1711,1712) with mi > 0, or with 777,2 > if m\ = 0. In the 
continuum limit of the single layer model, this assignment of the weight is possible by 
a certain modification of the partition function [12], which will be used in Section [3j In 
view of constructing the partition function basis (see the first footnote for single-layer), 
it seems to be necessary to include this type of distinct weight n in the definition. 

Here we give a general consideration on the distinction of the physics according to 
the value of A. Consider an arbitrary loop configuration on the layer A, then the weight 
for adding new non-contractible loop on the layer B will be, depending on the possible 
orientations, w + = n + X E ° or W- = n^\~ E °, where Eq G Z is the additional intersection 
number. If |A| = 1, we have w + = w*_, and both orientation of the new loop occur in 
amplitudes with the same modulus. On the other hand, if |A| 7^ 1, this symmetry is 
violated; one of the two possible orientation is preferred to the other. In this paper, we 
only consider the "symmetric" case: |A| = 1. Thus, we set 



for r G R. The scale of the parameter V is chosen for the later convenience, and the 
rescaled coupling g will be defined in (1231) . 

2.2. Counting the number of intersections between non-contractible loops 

Our purpose in introducing the signs to the intersections is to set the weight induced by 
the interlayer-coupling independent on a local fluctuation of loops. This is guaranteed 
by the invariance of the intersection number E under 




(8) 



oriented 
loops 



A = e 2 ^ r , 



(9) 
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• (i) an addition of contractible loops (Fig. E2-(a)), 

• (ii) a local deformation of loops (Fig. H3-(b)). 

Actually, it is these properties that enable us to decompose the bosonic field into local 
and global part and to do a path integral calculations in Section [3l 




Figure 3. (a) An addition of one contractible loop, (b) A local deformation of a loop, 
(c) Before the deformations. The opposite sides of the rectangle are identified. Since 
intersections indicated by the black and white dots have relatively opposite signs, the 
total number of intersections E is common (equal to five) for all the three cases. 



Now let us count the total intersection numbers E between loops belonging to the 
distinct layers. Firstly, one may forget about contractible loops by the invariance under 
(i). The result can be stated in terms of the pair of total winding numbers (l A , -m^) 
and (Jb, —ttlb) due to the presence of the several non-contractible loops on each layer. 
Below, we show the following: 

E(l A , -m A ; h, -m B ) = det ^ ™ A ™ B \ . (10) 

Since loops on a honeycomb lattice should be non-self-intersecting, only one type of 
homotopy class is allowed on each layer (otherwise they intersect within the same layer). 
This reduces the problem to the calculation of the intersection number between more 
fundamental objects. Namely, they are the torus knots characterized respectively by 
the winding numbers (l' A , —m' A ) and (l' B , —m' B ) which are coprime 

m' A M' A = m' B M' B = \, (11) 

where m A I stands for the greatest common divisor (GCD) of m and I; one should 
merely multiply ah to E(l' A , —m A ; l' B , —m' B ), if the configuration consists of the total 
winding number with m A Al A = a and itlb Al B = b. This procedure is compatible with 
the bi-linearity of the determinant. Thus, we can assume (fTTj) without loss of generality. 

To do the rest of the calculation, we consider an infinite plane with a coordinate 
(£, rj) in which the points differing by integer coordinates are identified and use it to 
parametrize the knots. Each of them is then represented by a "closed" curve in the 
unit square (£,77) G U = [0, 1] x [0, 1]. By the invariance under (ii), we may stretch the 
curve to form a set of pararell lines; these are periodic geodesies on the torus. We locate 

% The reason for the minus before m will become clear when we define the frustration in Section 13.11 
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these to pass through the origin. Then the locations of the intersections in U and their 
mirrors outside U are given by the solutions (£, rj) of 

m' A {Z-Z ) + l , A {7)-7] )=0 (12) 
rri B (Z-£ Q ) + r B (r}-Tj Q ) = 0, (13) 
where each set of (£o,r} Q ) G Z 2 and (£o>?7o) G Z 2 defines two lines and their unique 
intersection point. Since, by the Bezout's lemma, (II lj) implies that there exists a pair 
of integers (£o,r} ) G Z 2 such that m^£o + ^a^o = a f° r an y integer a and so does 
{Co,Vo) G Z 2 such that m' B £ + /^r/ = P for any integer (3, one can equivalently count 
the solutions of, 



m. 



m 




G Z 2 



(14) 



M = r 

The set of the solutions (£, 77) forms a lattice, which can be linearly mapped from the 
unit square lattice Z 2 by the multiplication of the inverse of the matrix in ( Ti~4l) . By 
taking the Jacobian factor into account, one can count the number of the solutions in 
U . We have thus shown the formula (jTOjl . 





Figure 4. (a) The streched configuration that is topologically connected with the 
configuration in Fig. |3]-(c). Again the number of intersections is five; the intersection 
at the top-right and bottom-left shoud be counted as one-half for each because 
of the identification. The winding numbers are given by {—IajTUa) = (3,1) and 
(— /s,ms) = (1,2). (b) The fundumental region U and its mirror on (£, 7y)-plane. 



3. Continuum partition function 

We first recall the Coulomb gas partition function of the 0(n) model [12J. Using the 
determinant formula ( flOl) . we formulate the continuum partition function of the layered 
model (|HJ). Here the moduli parameter of the torus is denoted by r = u 2 /uji = x + iy 
where uj\ and w 2 are the two basic periods. 

3.1. Single layer 

We introduce a compactified bosonic field (p(z) governed by the action 

A[<p] = ± J d 2 z (Vv^)) 2 . (15) 
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This field can be considered as a continuum version of the height function h(ri); the 
latter is defined on the dual lattice {rj} and changes its value by <ft = ±2irf when 
goes across the loop. The sign is again dependent on the orientation of the crossing; 
it is determined in a way parallel to the definition of op below dSJ). Since loops are 
closed, the correspondence between the loop and height configuration is unique up to a 
constant shift of the latter. On the torus, a configuration of the field f may belong to 
non-trivial topological sector called the "frustrations" along the periods (ui, 0J2), which 
are denoted here by (5i<f, 62^) = 2iif(m, I) with (m, I) e Z 2 . Now, in the continuum 
limit, the partition function j5| on a torus becomes 

ZgaussO, f) = Yl Z ™Mi f), (16) 



iii.t 



Z (g)= / V<p e~ A ^ = J ^\ V (q)\-\ (18) 

y 



Z m A9, f) = [ V<P e" AM = Z (gf) exp (- ^P\mr 

J(S 1 <p,8 2 ip)=2TTf(m,l) V V J 

where the field if is divided into global and local part: if = <f mt i + <fo- Quite importantly, 
the local fluctuation <f can be integrated out, giving Z (g) as [9] 

-AM = Jl\ n („\\-2 

where we have used the variable q = e 2mT and the Dedekind eta function 

00 

ri(q) = q* l[(l-q n ). (19) 

n=l 

In the standard convention, <pQ = it (f = 1/2); the coupling g is determined by El 

n = —2 cos ng. (20) 

Such description by the local action (|T5|) corresponds to the vertex rule <Q. Thus, the 
weight for the non-contractible loop (without orientation) is automatically h = 2. A 
novel observation [12] is that one can modify the weight for non-contractible loop as 

n = 2 — > n = 2 cos neo (21) 

by exploiting the identity: 

(2cosvre ) M = ]J ( ^ e^ ] = ]Tcos ( vre f>, ), (22) 

3=1 \e =±l j {e 3 } \ j=l J 

where {ej} represents 2 M ways of the positive or the negative orientation of M non- 
oriented non-contractible loops. A presence of the loop which is characterized by the 
coprime winding numbers (l',—m') leads to the frustration (5\<f, 62^) = 2nf{m',l'). 
Thus, the net number of oriented loops is related to the frustration as 

M 

£' 

i=i 



mAl. (23) 



+ The value of g is fixed from the requirement that the screening field e ltp /f, which originates from the 
discreteness of the height function h(ri) £ 2nfZ, be marginal [13] . This leads to the relation g = 1 ±x 
with x defined in 
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As the result, the partition function is given by 

Z o(n) = ^2 Zm > 1 cos ( 7re o m A > ( 24 ) 

m,l 

where we have introduced, for convenience, the rescaled coupling 

9 = f , (25) 

and have used short-hand notation Z m j (g) = Z m ^\ (g, 1) and the relation Z m j (gf 2 , 1) = 
Z m ,i (9 J) with/ = 1/2. 

3.2. Bi-layer 

From the lattice expression ([8]), definition (Q, and formula (fTUj) . the partition function 
of the layered 0(n) model can be written as 

Z = ^2 Zm A ,i A (g) cos (ne m A A l A ) C{m A , l A , m B , l B ; T) 

m A ,l A m B ,l B 

■ Z mB ,i B (g) cos (ne m B A l B ) (26) 
where the inter-layer coupling is given by 



C(m A , l A , m B , h] r) = exp 



27cigT ■ det 



m A m B 



(27) 



At the level of the action, this term would correspond to the coupling 

SabM = ! d 2 zh- (V VA (z) x Vip B {z)) , (28) 

In J 

where h is the positive unit vector perpendicular to the ^-complex plane. Since the 
integrand has the canonical dimension two, a possibility is that the presence of the 
inter-layer coupling leads to a marginal deformation of the decoupled model (r = 0). 
In general, an invariance of partition functions under the modular transformation: 

T ^£I+& fa M e5L(2 , Z) , (29) 
cr + d \ c a I 

imposes strong constraints on the operator content [8]. In this respect, our expression 
(126"]) of the bi-layer system automatically satisfies the modular invariance. To see this, 
first recall that any modular transformation (129 j) can be expressed as a finite composition 
of the two basic transformations T : r — > t + 1 and S : r — > — 1/r. The modular 
invariance of ( fl6l) follows from the transformation properties of the frustrations (g- 
dependence being suppressed): 

T : Z m j (r + 1) = Z m ^ m (r) , (30) 
S: Z m .i f-^j =Z^ rn (r), (31) 

* The author thanks J. Cardy and J.L. Jacobsen for comments concerning this point. The author also 
notes that a similar (but different) action was used to study the conformal theories which posses W n 
symmetry 1 . 
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and from the basic idea of the path-integral that the sum should be taken over all the 
possible configurations (here, taken over all the frustrations). Since the GCD remains 
invariant under these transformation: 

m M — m A {I — m) — (-/) A m, (32) 

the partition function of the 0(n) model (12^1) also has the modular invariance. Further, 
the determinant has the analogous invariancej^] which may be represented as 

M xL = M x(L- M) = (-L) x M (33) 

if the determinant in f )27l) is regarded as the third component of M x L with M = 
(mAytTiByO) and L = (Ia,Ib,®)- Thus, under the transformations, each term in the 
sum ( 126]) changes into another term with the common extra factor, compared with ffl6|) , 
determined only from the GCD and determinant of the frustrations. This establishes 
the modular invariance of the partition function of the layered model f l26|) . 

In order to see the operator content, we need to Fourier transform the sum ( 1261) into 
a series in q and q. The similar transformation for the partition function ( TT6|) is easily 
done by the Poisson ressumation. Intuitively, this task becomes rather non-trivial for 
the sum with the GCD, (12^1) and ( )26l) . because the GCD function mAl has no apparent 
periodic structure on the (m, /)-plane; its distribution has rather recursive (self-similar) 
and irregular pattern (being recursive because all the prime numbers play equivalent 
and relatively independent roles in determining the GCD, and irregular because they 
appear irregularly). It is possible, however, to take advantage of the recursive nature by 
allowing deliberate overlaps and to express the sum ( 12^1) as a superposition of sums over 
the periodic structures (m, I) e kl? . In order to do this, we use the following formula: 

^2 Z m,i (g) cos (ire m A /) = ^ ^2 ^ ) cos ( 7re o rf ) /2 Zm > 1 ( 34 ) 

m,l d>0 k: d\k ^ ' k\(m,l) 

where the Mobius function is defined by 

{1 if n = 1, 

(— l) k if n is a square-free integer with k distinct prime factors, (35) 
if n is not square-free. 

Employing such technique based on a version of the Mobius inversion formula |40j, 
the partition function of the 0{n) model was evaluated [39]. Changing the order of 
summation over d and k in and then using that formula twice, we obtain 

fcyt>0, k B >0 d,A\kA, d B \k B 



y~] ^ Z mAt i A C(mA,lA,m B jB;r)Z mBt i B . (36) 

kA\(m A ,lA) k B \(m B ,l B ) 

jj The invarince of the determinant under the general modular transformation can be also seen by the 
fact that the frustration vectors (mA,lA) and (jubtIb) transform by multiplication of the identical 
SX(2,Z) matrix, which has the determinant equal to one. 
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Note that the sum in the second line factorizes as below: 

Zm Ai i A C(m A , I a, m B , h; ^)Z mB ,i B 

kA\(m A ,lA) k B \(m B ,l B ) 

= ^ e -^[( m A x + l A) 2 +rn\y 2 ] e 2TTigT(m A l B -l AmB ) e ~^[{m B x+^ 

I A- ^a\Ia, 
l B '- k B \l B 

= s£\-m B )Sl B B \+m A ), 
where the last two factors are given by 

S^\- mB ) = ^2 e~^ [(U+m A x) 2 +m A y 2 ]+2nigri A -(-m B ) 
I A- Wa 

l B : k\l B 

Then we shall perform the Poisson resummation to each quantity. It is sufficient to see 
this for, say S^\— m B ); this goes as: 

5<( A )(_ mB ) — ^2 e -^[( l A+m A x+irm B y) 2 -2irm A m B xy+(r mB y) 2 +m 2 A y 2 ] 



(37) 

(38) 
(39) 



l A : k\l A 



= J—y 




-wy 



-ny 



^ I +g{Ym B ) 2 +gm 2 A 



./ m AX-\-iTm T^y \ 



m A m B x 



^+g(rm B ) 2 +gm 2 A+ 2 r -^ 



Ilk 



_^2nix(^ —j^--\-gTrn, A rn B ) 



(40) 



Comparing this with the expression q h q h = exp [— 2ny(h + h) + 2mx{h — h)\ , we obtain 



Q^(—rn n \ - n Hk,p,m A ,m B )-h(k,p,m A ,m B ) 
D fc I m B) - \ p lh 2 q ' 



with 



h(k,p, m A ,m B ) 



h(k,p, m A ,m B ) 



1 / p 



4 \^k 
1 / p 



g (m A + Tm B ) 



- V9 ( m A - rm B ; 



(41) 



(42) 



(43) 



4 \^gk 

The magnetic charge of the layer A and that of B are mixed by the presence of the 
interlayer coupling V . If we define 

A(p A /k A ,p B /k B ,m A ,m B ) = h(p A /k A ,m A ,m B ) + h(p B /k B ,m B , -m A ), (44) 
A(p A /k A ,p B /k B ,m A ,m B ) = h(p A /k A ,m A ,m B ) + h(p B /k B ,m B , -m A ), (45) 

then the scaling dimension and conformal spin yield 



A + A = 



1 



£(if + j)+«(i + r 2 ) K + •»«] + ( 



PA PS 
— m B - m A — 
k A k B 



A- A- Va1Ua + P b7Ub 



k, 



k R 



(46) 
(47) 
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where the arguments of A and A are the same as the one written in 
Substituting equations f|T8|) . f )37|) . and fj4Tl) into fl36l) we obtain 



and (US]). 



E E 



E E «V. 



cos (7ie dA) cos (7ie dB) 



(4* 



vriA'- kA\m,A,PA&'i 
1TIB'- kB\lTlB PB& 1 



In order to see the operator content, we should evaluate the coefficient of q A q A with 
fixed sets of values (A, A). Since these values depend on the electric charges — Pi/ki 
{i = A, B) which are invariant under (pi, ki) — » (cpi, cki) for positive integer c, we would 
like to group these terms, which gives the same values of (A, A). 

For this purpose, it is natural to decompose the sum (148]) into terms with rriArtiB ^ 
and mxm B = 0, which are respectively denoted by Z^ and Z so that 

\ V (q)\ 4 Z = Z^ + Z 



(49) 



For the sum Z^q , we apply the formula derived in |Appendix A 



k>0 d\k 



cos(vre rf) ^ E^(? 



te- 



rn 



m: k\m,m^0 p& 

= EE A ^ E f{v m . 

m^O k: k\m p: pAfc=l 

where / [z,Tnj is an arbitrary function (thus, not necessarily even in the variable 
the notation m = \m\ is used for readability Here, A(k,m) is the multiplicity in the 
single layer 0(n) model given by the finite sum: 



(50) 



m 



A (k, m) 



m 



E c * 

d: d\rh 

where C a (b) is the Ramanujan sum defined by 



cos (7re c?) 



Caib) = 



(51) 



(52) 



hf\a=l 



Then we have 



Z 



q^q \ k l ,k 2 , L ' z ) 



(53) 



n E E E 

1=1,2 rriij^O kii ki\rhi p%: PiAki=l 

Next, we consider the sum ^ which consists of the terms with which at least one of 
■itla and tub vanishes. Taking the terms with mg = and itla ^ 0, for example, we 
see from ( |26|) that the sum over I a reduces to the nonzero-magnetic-charge sector of the 
single layer 0{n) model, which can be evaluated using the formula (134"]) . Further, since 
the GCD is to be evaluated as A Ib = Ib, the sum over Ib is just 



EE 

s=±l l B 



-TY 

4(7 ^ ^ 

y s=±l Pb 



(54) 



Coulomb Gas Partition Function of a Layered Loop Model 



14 



Therefore we see that setting mg = amounts to a restriction on the electoric charge in 
( 146]) and (14"7]) to e B = p B /k B with k B = 1 and a shift in p B — )■ p B ± ^ + grm^. Taking 
the terms with tua = and the symmetry A (k, m) = A (k, —m) into account, we obtain 

P1,P2 

mj^O k: k\\m\ py. piAfc=l P2 

To summarize, our result on the partition function of the bi-layer model is given by (149]) 
with (1551) and (153]). 

4. Spectrum of the bi-layer models 

The first rather trivial observation is the following. When the inter-layer coupling V 
vanishes, the scaling dimension (T46|) becomes just a sum of two dimensions known in 
usual CG. Hence, the bi-layer model ( 148]) decouples into two independent 0{n) models: 

Z(g, e , T = 0) = [Z (n){g, e )] 2 (56) 
with Z O ( n )(g,e ) given by [H] 

\v(q)\ 2 Z o{n) (g,e )= ^ g A-A( Pl+? ,o,o,o) 

pi 

+ £ £ A(*,m) £ g V^ Am '°). (57) 

m^O k: k\\m\ pv p±Ak=l 

as it should be. 

From (jUJ), one can in principle read off the non-trivial selection rules in the bi- 
layer system. In particular, it is always possible to calculate the lower spectrum by this 
explicit formula. It seems, however, difficult for general cases to state the selection rules 
in simpler form e.g. in terms of few quantum numbers. We do not pursue this direction, 
and just present some examples in order to convey the flavor of how the selections occur. 

4-1- Spectrum in the bi-layered dilute/dense 0(1) models 

The single-layer 0(1) model which is realized at g — 4g — 4/3, eo = 1/3 [12] is the 
Ising model, which has the central charge c = 1/2. Hence, the bilayer model without 
interaction has c = 1. A partition function of CFT, in general, behave as 

Z ~ (qq)~^ . (58) 

By now, we use a combination which is suitable for analyzing c = 1 partition functions: 

Z cl (g,e ;T) = (qq)-^\r ] (q)\ A Z, (59) 

which would start from the terms of order 1 as the lowest-order term. In general, 
if a conformal field theory has a primary fields of dimension h, the theory also has a 
descendant field of the dimension h + n which is, in the algebraic setting [6], obtained by 
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applying the Virasoro generators L_ ki (hi ^ 1) to the primary such that n = k 1 + - ■ -+k m . 
As is well known, this property is closely related to the fact that the inverse of the eta 
function ( fl9l) is almost the generating function of the partition numbers d(n): 

oo oo 

^(g)- 1 = g -H JJ (1 - q n y X = q-^ Yl ( 60 ) 

n=l n=0 

Correspondingly, we should note that the coefficients of q h q h depends on the power of 
the eta function which we factored out. However, here we focus on the lower spectrum in 
( 15 9p and do not concern the higher order terms corresponding to the descendants. From 
now on, we study the coefficients which appear in fl59|) and call them "multiplicities" . 
Using g9]) with (J53J) and (J55D, we obtain 



Zcl Q>^°) = (l-2g + g 2 ) (l-2q + q 2 )+qq + 2^ _ giff) (qls - q\ 

+ 2 (qq)& + (gg)* + 2 (q? - g§J {(p - g§J + O (g 2 , q 2 ) , (61) 

where the symbol O (q 2 , q 2 ) represents the terms of higher-order than q 2 or q 2 . The result 
( |6T|) can be checked by using either the Ising partition function (IB.lj) or the orbifold one 
( 1B.6j) in Appendix B. Reflecting the non- interacting nature, the spectrum is given by 



the trivial sum of two elements chosen from the single Ising spectrum {0, 1/16, 1/2}. 

Next we turn to the interacting model. A simple, and yet non-trivial spectrum is 
realized at A = — 1 (r = 3/2). It reads, 

(1 X 3 \ 1 9 17 

-,-;-J = (1-g + g 2 ) (l - g + g 2 ) + 2 (qq) ™ + 2 (qq) * + 2 (qq) ™ 

- 2 - 2 (qq)T + 4 (q? - qi\ (q* - gi) + O (g 2 , g 2 ) . (62) 

Comparing the interacting model ( |62]) with the decoupled one fl6T|) . we notice the 
presence of the new characteristic dimension 3/8. This dimension also appears, for 
example, in the spectrum of the level-2 coset model [27J [281 122] 

s£(2) 2 ®s£(2) 2 
s£(2) 4 

or in the Ashkin- Teller quantum chain at special points [301 Ell 132] , where the spectrum 
can be related to the irreducible representations of the c = 1, M = 2 superconformal 
algebra. The dimensions which appear in the table in the coset model are 



^ 1 1 1 3 9 2 3". 



24' 16' 6' 8' 16' 3' ' 2 

Further, we remark that the three diagonal exponents contained in f l64j) but not in 
namely 1/24, 1/6, and 2/3 do appear also in the interacting model with g = 1/6, 



X In the c = 1, Af = 2 superconformal model, the two exponents 1/24 and 1/6 belong to the Ramond 
sector and to the Neveu-Schwarz sector, respectively [33], 
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e = 1/3 (n = 1), and T = 3 (A = exp [27ri^r] = —1), whose partition function reads, 
Zd Q'^ 3 ) = (gg)^(3-4g-4g + 6gg) -2 [(<?# + (<?#] (1 - g - <?) 

+ 0(g 2 ,g 2 ). 

The value g = 4g = 2/3 also leads to n = 1, which is the same as in (162]) , but belongs to 
a different branch in f[2"0l ; this value is known to correspond to the dense 0(1) model. 
The corresponding non-interacting model has the simpler function 

=4(g#(l-2g-2g + 3gg) + 0(g 2 ,g 2 ). (66) 

Here we just substitute the value for the dense model, although our argument (especially 
that in Section 12. 2j) is more reasonable for the dilute model. Note also that the lowest- 
order term in ( 165]) which is diagonal in q and q is (qq)* 1 , which means in our notation 
( 159]) that ( 165|) is related to a model with the central charge c = rather than to that with 
c = 1 if it is considered as a full partition function. It remains, however, the possibility 
that fl65|) forms a part of some other partition function with the central charge c = 1. 

In summary, the partition functions of the bi-layered dilute 0(1) models and the 
bi-layered dense 0(1) models both at the intersection weight A = — 1 are represented 
respectively by (|62|) and ( )65|) . and together reproduce all of the spectrum (l64"j) of the 
coset model ( 1631 . It would be gratifying if one could find more precise relation between 
the partition function of the coset model and special cases (which may include some of 
the examples in this section) of our modular invariant partition function. 



4-2. Flow of the spectrum in the bi-layered dilute 0(1) model 

The curves of conformal dimensions belonging to the sector Z^ and Z connecting the 
dilute 0(1) model at A = 1 (T = 0) and that at A = -1 (T = 3/2) are plotted in Fig. El 
and Fig. [61 respectively. As one can see from (jUJ), the spin A — A is independent of T. 
Accordingly, we choose to show the dimensions of spin zero (diagonal) operators. We 
plot A' = A — 1/24 instead of A itself. One can reproduce the diagonal terms in (I6T!) 
and that in (162|) . including degeneracies, by looking at the intersections of the curves to 
the left and to the right axis of both graphs. 

All the curves in Fig. [5] are parabolas, as is clear from (j4"B"]h They are folded due 
to the symmetries of the partition function that follow from the definition ([9]): 

r ->■ -r, r -> r + 3. (67) 

In addition to these, as in Fig. [6j we have an extra symmetry 

r^r + | (68) 

in the sectors Z . Formally, this can be understood by the presence of the half-divided 
charge eo/2 in f )55|) . The four flat lines (A' = 0,1/2,1, and 3/2) correspond to the 
operators with vanishing magnetic charges and mg. As a result of (168]) . all the 
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Figure 5. Flow of the conformal dimensions A' = A — 1/24 as functions of T in the 
diagonal sectors of Z^ in the model with g = 1/3 and 6q = 1/3. Each curve carries a 
definite multiplicity (solid: +1, thick: +2, dashed: -2). 




0.0 0.3 0.6 0.9 1.2 1.5 



r 



Figure 6. Flow of the conformal dimensions as functions of T in the diagonal sectors 
of Z in the same model (g = 1/3 and eo = 1/3). Again each curve carries a definite 
multiplicity (solid: +1, thick: +2, dot-dashed: +3). 



difference of the operator content of the model at A = 1 to that of the model at A = — 1 
are contained in the sector Z^q (Fig. [5]). 
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Thus, let us concentrate on the flow in the sector Z^ . By recognizing that these 
parabolas occur in sequence, we may state an interpretation of this flow as follows. If 
we look at the graph up-side down, these curves are regarded as some trajectories of 
balls in a space surrounded by the walls at T = and at T = 3/2 in the presence of a 
uniform gravitational field. For example, if we are at (r, A') = (0, 1/8) and throw a ball 
in the horizontal direction to the wall at T = 3/2, then it falls down to A' — > oo as it 
reflects alternatively at T = and at T = 3/2. If we tune the initial velocity, it goes as 
(0, 1/8) — > (3/2, 1/2) — > (0, 13/8) which in fact appears in the graph. Another 

trajectory is (3/2,3/8) — > (0,21/16) Every collisions are perfectly elastic due to 

the symmetry (1671) . Hence, we have a law of the reflection in the spectrum. 

Each falling trajectory carries a certain positive or negative multiplicity and implies 
the possibility of the operators of higher scaling dimensions at any T. For fixed T, it 
defines a sequence of "descendants"^. At T = 0, for instance, the descendants of the 
field of dimension A' = 1/8 may have dimensions {13/8, 49/8, 109/8 •• •}. Finally, we 
remark that cancellations of multiplicities may happen at particular values of T where 
multiple trajectories cross each other. As we can see from (I6T1) . there are no diagonal 
terms with A' = A' = 13/8, or 21/16. Complete cancellations in fact occur both at 
(0, 13/8) and at (0,21/16). This illustrates a non-trivial aspect of the selection rules. 

5. Conclusion and Discussions 

We have introduced a bi-layered loop model in which the inter-layer coupling depends 
on the signed number of intersections between loops living on the different layers. The 
model is first defined microscopically on the lattice with the periodic boundary condition, 
and then the corresponding continuum partition function is proposed based on the 
idea of the Coulomb gas. Within our definition of the inter-layer coupling, the local 
fluctuation decouples from the topological one. Accordingly, we are left with the sum 
over the topological sectors of two compact bosonic fields. 

The partition function has the modular invariance. After the Poisson resummation, 
we have the explicit expression which reveals the operator content. Each coefficient 
involves a combination of the Ramanujan sums, and contains non-trivial selection rules. 
We analyze few examples of the curves of scaling dimensions induced by the inter- 
layer coupling (Fig. [5] and [6]). The layered model continuously changes itself from 
the decoupled model (A = 1) into the non-trivial coupled model (A = —1), where the 
operator content is reorganized and the model has the simple spectrum f l62l . 

In the definition of our layered model on the lattice, loops from different layers do 
not overlap but instead intersect each other. This feature is unique to our model. To 
compare, we mention the coupled loop model whose self-dual critical point was studied 
by algebraic methods [5] and the models that arise from the replica approach to the 
disordered models such as the random-bond g-state Potts model [U [2] or the disordered 

§ Here, in general, difference of the dimensions are restricted to neither integers nor half-integers; they 
are not the usual descendants well-known in the conformal, or superconformal theories. 
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0(n) loop model [3J. In the former model, two Potts models, described respectively by 
two independent Temperley-Lieb algebras, are considered on the same square lattice and 
coupled via the weight for the local overlap of loops. Similarly, for the latter models, the 
averaging procedures lead to the models of the replicated M-layers which are coupled 
through the overlap of loop segments (energy operators), that is, renormalized version 
of the bonds Js ■ s where s is the microscopic spin. Although, at this stage, we could 
not relate our model and these other coupled loop models, we briefly discuss a possible 
approach to the latter, disordered models by using the method in this paper. 

In the replica method, the intermediate layered models were analyzed for general 
M G N by conformal perturbation theories, and then M — >■ limits should be taken. 
In this way, the two lines of the non-trivial fixed point emanating from the Ising point 
q = 2 and n = 1 are found |H El E]. Along these critical lines, the effective central 
charge and, thermal and magnetic eigenvalues are estimated. They are also numerically 
checked, for the random Potts case, using the transfer matrix [35j EH]. Now our model 
can be generalized to the one consisting of M-layers. Since it would be again described 
by a quadratic action similar to that used in this paper, it seems to be straightforward 
to obtain the partition function of the model for any M G N. It would then give some 
insight on the replica limit if the analytic continuation in the variable M and studying 
its M — > limit could be possible at the level of the exact partition function. 

Further important extension is within the bi-layered model but is a generalization 
including several number of charges (i = 1,---,N), or equivalently, with several 
weights hi for non-contractible loops. Other than at h G N, this is necessary for 
obtaining integer coefficients. In the single-layer case, {hi} is taken to be the set of the 
eigenvalues [17] of the adjacency matrices associated with the A-D-E Dynkin diagrams. 
In this way, it would be possible to examine whether we have a discrete series of the 
models similar to the unitary minimal model [7] in the single-layer case or not. 

Finally, we emphasize that the virtue of the CG method is its directness. Although 
we have only dealt with a specific case in multilayered system, the continuum partition 
function was derived from the lattice formulation by a straightforward path-integral. 
We hope that our slight extension in two dimensions from single layer to multilayer will 
stimulate further progress in this direction. 
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Appendix A. Derivation of the formula ( 1501) 



We here present the derivation of the key formula ( 1501) for completeness although a 
similar calculation has been given in the literature [39]. Let F be the left hand side of 
( l50l) . Changing the order of the summations and then rewriting by the reduced fraction 
P/K yields (below rh denotes \m\ for readability), 

^EE^w E £/(H ( a -d 

fc>0 d\k m: k\m p& 

m^O 

^EEEx^wE/l-) < A2 > 

m^O k: k\rh d\k pgZ 

= EE E Ef™w E /(£.-»)■ < A3 > 

More convenient expression is obtained by setting k = rh/t: 

F = E E E E^m E /(£.»•) < A4 ) 

m^OK: K\mp. K\fd: d\f ' P: PAK=1 V ' 

= E E E E ^f^ C03(wi) £ /(£, m ). (a.5) 

In order to obtain ( 1A.5I) . we have changed the order of the summation over t and that 
over d; the region of the summation is visualized in Fig. IA1I After evaluating the sum 
over t explicitly by applying the formula [40] for the Ramanujan sum: 

v(^)t = C a (b), (A.6) 

t: t\(aAb) 

we obtain the result in (I5U1) . 




Figure Al. The change of the order of the summations. The summation is over the 
trapezoid in both figures; a log-log scale is adopted so that the relation td = rh holds 
along the hypotenuse of the right triangle. 
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Appendix B. The Ising and orbifold partition functions 



The partition function of the critical Ising model (the dilute 0(1) model) is, 

4 

■^Ising — 2^ 



j=2 







J 

277 



(B.l) 



where 9j (j = 2, 3, 4) are the classical theta functions (at z = 0) defined by 



e 2 (z; t) = J2 q*^*? e 2n < n+ ^ z = e niz q l s f[ (1 - q n ) (l + e 2niz q n ) (l + e~ 2niz q n ) , 

nGZ n=l 
oo 

q^ e 2ninz = Y[ (1 - <f ) (l + e 2niz q n '^A f\ + e " 27r V*) , ( B .2) 



n=l 



-2niz „n~ 



e A ( Z] r) = ^g^V^+s) = J[ (1 - ? n ) (l - e 2niz q n ~ ^ ) (1 - e 

n6Z n=l 

We need the square of ( IB. II) . Here, it may be instructive to recall the partition function 
derived from the orbiforld boson (or that of the critical Ashkin- Teller model [331 132]). 
which includes the decoupled two-Ising models as a special point. It is given by 

i 



Zorb(g) = ^Z geMBB (g, 1) + ^2 



3=2 



(B.3) 



Formally, one can see the last sum corresponds to the cross term in the square of fIB.ll) 
by noting the identity between the eta function (|T9j) and theta functions (1B.2j) : 

2 V 3 = 9 2 9 3 9 4 . (B.4) 

Physically, the sum in f ]B.3j) has the meaning of the sum over the triplet of the spin 
structures (they exchange each other under the modular transformations) [9] 



J AP 



04 



J PA 



J AA 



(B.5) 



where the subscripts corresponds to the boundary conditions along the two periods of 
the torus; A and P represent anti-periodic and periodic, respectively. At the special 
point g = 1/2, (IB.3j) indeed reduces to the square of (IB. II) : 



Z rb ( g — - 



v Ising> 



(B.6) 



as stated. 
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